Spectrum of weak magnetohydrodynamic turbulence 
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Turbulence of magnetohydrodynamic waves in nature and in the laboratory is generally cross- 
helical or non-balanced, in that the energies of Afven waves moving in opposite directions along 
the guide magnetic field are unequal. Based on high-resolution numerical simulations it is proposed 
that such turbulence spontaneously generates a condensate of the residual energy E v — Eb at small 
field-parallel wave numbers. As a result, the energy spectra of counter-propagating Alfven waves 
are not scale- invariant. In the limit of infinite Reynolds number, the universality is asymptotically 
restored at large wave numbers, and both spectra attain the scaling E(k) cx kj 2 . The generation of 
condensate is apparently related to the breakdown of mirror symmetry in non-balanced turbulence. 

PACS numbers: 52.35.Ra 



(N ; 

O: 

Oh; 

o 



> 
in 

o 

O 



Introduction. — Magnetohydrodynamic (MHD) turbu- 
lence naturally occurs in a variety of plasmas, ranging 
from the interstellar medium to the solar wind to labora- 
tory fusion devices. When the compressibility effects can 
be neglected, the MHD equations take especially simple 
form in the so-called Elsasser variables, 
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where z = v ± b, v is the fluctuating plasma velocity, 
b is the fluctuating magnetic field normalized by ^/AiFpo, 
va = Bo/V47TjOo is the contribution of the uniform mag- 
netic field Bo, P = (p/po + b 2 /2) includes the plasma 
pressure p and the magnetic pressure, po is the constant 
plasma density, represents the mechanisms driving 
turbulence, and small dissipation due to viscosity and re- 
sistivity is neglected. In the absence of dissipation, both 
energies E + = (|z + | 2 ) and E~ = (|z~| 2 ) are conserved, 
which is equivalent to the conservation of total energy 
and cross- helicity 

The linear terms, (va ■ V)z ± , describe advection of 
Alfven wave packets along the guide field, while the 
nonlinear interaction terms, (z^-V)z ± , are responsi- 
ble for energy redistribution over scales. Depending 
on the driving force, turbulence exhibits either weak or 
strong regime in a certain range of scales. Denote b\ the 
rms magnetic fluctuations at the field-perpendicular scale 
A cx l/k±, and assume that the typical field-parallel wave 
vector of such fluctuations is fey . Then the turbulence is 
weak when the linear terms dominate, kuVA 3> k±b\, and 
it is strong otherwise. 

Weak MHD turbulence may be present in laboratory 
devices and in the solar wind, as indicated by energy 
spectra somewhat steeper that the Kolmogorov one 
Moreover, for a general driving force weak turbulence 
regime precedes the development of strong turbulence, 
in a sense that as the scale of fluctuations decreases the 
turbulence eventually becomes strong. 

When the nonlinear interaction is absent, the solution 



of {T]) is an ensemble of shear- Alfven and pseudo- Alfven 
waves propagating along the guide field Bq with the ve- 
locities ±va- The small nonlinear terms then can be 
taken into account perturbatively, and the spectrum of 
turbulence can be derived using the general methods of 
the theory of weak turbulence; for a review see 0, 0]- 

Spectra of MHD turbulence were first studied by Irosh- 
nikov [3] and Kraichnan Those early works realized 
the role of the guide field in mediating the turbulent 
cascade, however, they assumed the small scale fluctu- 
ations to be isotropic. Over the years, the assumption of 
anisotropy proved to be incorrect [e.g.,[H, Anisotropic 
spectra of weak MHD turbulence were addressed by Ng 
& Bhattacharjee Q and Goldreich & Sridhar Q based 
on dimensional arguments, and a comprehensive analytic 
framework was developed by Galtier et al [9(. The lat- 
ter theory derives the kinetic equations for evolution of 
the spectral energies e ± (k) = (|z ± (k)| 2 ), and has the 
following main results. 

First, the spectral energies are transferred in the direc- 
tion of large k± , and the universal regime of weak turbu- 
lence is established at k± ^> k\\ ■ In the universal regime, 
the dynamics of shear-Alfven modes decouple from the 
dynamics of pseudo- Alfven modes, which are passively 
advected by the shear-Alfven ones. In what follows we 
shall consider only shear-Alfven waves and keep the same 
notation e ± (k) for their energies. It is also customary to 
use the phase-space- volume compensated spectra defined 
as ^(k) cx e ± (k)fc J _. 

Second, the predicted spectra of weak MHD turbu- 
lence are not unique, but form a one-parameter family, 
E ± (k±) cx k^ 2±a , with -1 < a < 1. The solutions 
with a / correspond to unequal fluxes of the i? ± (fc)- 
energies over scales; we denote these fluxes e + and e~. It 
can be further demonstrated that in this case the ampli- 
tudes of the energy spectra are different, e.g., E + > E~ 
if e + > e - . The spectral parameter a is uniquely defined 
once the ratio e + /e~ is specified. In the balanced case 
the energy spectrum is E ± (k±) cx kj 2 0,S, 
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MHD turbulence in nature and in the laboratory is of- 
ten driven by localized sources (solar wind, antennae 
localized instabilities, etc.) and therefore it is gener- 
ally non-balanced. There is however a more fundamen- 
tal reason for imbalance in MHD turbulence. Numeri- 
cal simulations demonstrate that even when MHD tur- 
bulence is balanced overall, it gets spontaneously non- 
balanced locally, that is, it consists of domains of pos- 



ii 



This do- 



itive and negative imbalance 
main structure is apparently related to the conserva- 
tion of cross-helicity in ideal MHD, and non-balanced 
turbulence is also called "cross-helical." Non-balanced 
MHD turbulence has recently attracted considerable in- 
terest 0, Ei EI EI EE E3 . 

In the present work we study non-balanced MHD tur- 
bulence in a series of high resolution numerical simula- 
tions. The results reveal puzzling contradictions with the 
theory. While in the balanced case the numerics confirm 
the analytic prediction E(k±) oc kj 2 [cf 0, EH, in a 
general non-balanced case simulations disagree with the 
theory. When we drive turbulence with unequal rates e^, 
the resulting spectra E ± turn out to be not defined by 
the ratio e + /e~. Rather, they depend on the Reynolds 
number and approach fcj 2 at large kj_ as the Reynolds 
number increases. 

To resolve this contradiction, wc propose that driven 
weak MHD turbulence generates the residual energy 
condensate (z + (k) • z~(— k)) = v 2 — b 2 ^ at fen = 0. 
This condensate has been assumed to be zero in the 
standard derivation 0, [H, Q, therefore, its presence in 
our model requires an explanation. The Alfven wave 
fluctuations obey v = ±b, in which case the residual 
energy vanishes. However, at fen = fluctuations are not 
waves, and the Alfvenic condition should not be neces- 
sarily satisfied. We further propose that the generation 
of the condensate is a consequence of the breakdown of 
the mirror symmetry in non-balanced turbulence. 

Kinetic equations for weak MHD turbulence. — In this 
section we discuss the predictions of the standard model 
for weak MHD turbulence developed in @, [1, @, El ■ The 
kinetic equations for the evolution of the shear- Alfven 
energies, derived by Galtier et al [9(, have the form: 
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: (q) [e ± (p)-e ± (k)] 5( gil )d km ,(2) 



where we use the short-hand notation 
{■k/v a ){\l± x q_L) 2 (k_L ■ p 1 _) 2 /(k 2 J _p 2 L q 2 1 _) and d ktPq 
<5(k p — q)d 3 pd 3 q. The derivation assumes that in the 
zeroth-order approximation, only the correlation func- 
tions e ± (k) are non-zero. As shown in the system (|2|) 
has a degeneracy: the right hand side integrals vanish for 
any solutions of the form 



with arbitrary functions g^(k^) and — 1 < a < 1. The 
degeneracy is removed by matching these solutions with 
the boundary conditions, that is, forcing and dissipation, 
[e.g., E3|- To match with the forcing, one notes that 
different energy spectra correspond to different energy 
fluxes, supplied by the large-scale forcing, such that 
a is uniquely found if the ratio e + /e _ is specified. One 
can show that the solution with the steeper spectrum 
corresponds to the larger energy flux . 

The large-scale boundary conditions fix the slopes of 
the energy spectra, but do not fix their amplitudes. To 
fully remove the degeneracy, one further argues that at 
the dissipation scale the balance should be restored, that 
is, e+(fc) should converge to e~(fc). This "pinning" effect 
was first pointed out in 1 131 and it physics was discussed 
in greater detail in (9Lll4Lll5|. The pinning effect is indeed 
observed in our simulations presented below. 

According to the above picture, if the rates of energy 
supply are fixed, then the slopes of the energy spectra 
e ± (fc) are fixed as well. If the dissipation scale is now 
changed, the amplitudes of the spectra should change 
as to maintain the specified slopes, and to make them 
converge at the dissipation scale. This conclusion, 
although consistent with equations ([2]), seems to be 
at odds with the common intuition about turbulent 
systems, which suggests that small-scale dissipation 
should not significantly affect the large-scale fields 
subject to the same large-scale driving. This seeming 
contradiction motivated our interest in the problem. 

The numerical method and the results. — The universal 
properties of MHD turbulence with a strong guide field 
are accurately described by neglecting the field-parallel 
components of the fluctuating fields, associated with the 
pseudo- Alfven mode [TEHIl- By setting zjj 1 = in equa- 
tion ([1]) we obtain the closed system of equations 



= -V,P 



-ff+j/V 2 z ± , (4) 



--(k) = g ± (k ll )k[ 



-3±a 



(3) 



in which dissipation terms have been added, and we as- 
sume that viscosity is equal to resistivity. This set of 
equations is known as the Reduced MHD model origi- 
nally developed for tokamak plasmas 2lj, [22J, and of- 
ten used in numerical simulations of different regimes of 
MHD turbulence. Depending on the spectral properties 
of the driving force, this system can describe either weak 
or strong MHD turbulence [ToL [20I] . 

We employ a fully dealiased Fourier pseudo-spectral 
method to solve equations (|4]) with a strong guide field 
{vA/vrms ~ 5) in a rectangular periodic box, with field- 
perpendicular cross section L\ = (27r) 2 and field-parallel 
box size Lu = 5L±. The choice of a rectangular box, as 
discussed in [10fl . allows for correct description of long- 
wavelength and low-frequency fluctuations. 
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FIG. 1: Upper panel: The spectra of balanced weak MHD tur- 
bulence. The solid line is E + (k±), the dashed line is E~ (k±); 
Re = 6000, resolution 1024 2 x 256 points. Lower panel: The 
spectra of non-balanced weak MHD turbulence, with the im- 
balance parameter 7 = 0.17. The solid lines denote E + (k±) 
and E~(k±), Re = 4500, resolution 1024 2 x 256 points. The 
dashed lines show the same fields for Re — 2000 and resolu- 
tion 512 2 x 256 points. The insert shows the corresponding 
energy dissipation rates. 



The z + and z waves are driven independently by 



Gaussian random forces 
F±^2 



and f - , with the variances 



a = ((f|_ ) 2 ). The imbalance is measured by the param- 
eter 7 = (er + — <t~ )/ (<t + + er~ ) . To ensure that the turbu- 
lence is weak, the forces have a broad fcii spectrum. They 
are applied in Fourier space at wave-numbers 1 < k± < 2 
and (2w/Ln) < kn < 16(27r/-Lii). The Fourier coefficients 
inside that range are independent Gaussian random num- 
bers with the amplitudes chosen so that the resulting 
rms velocity fluctuations are of order unity. The individ- 
ual random values are refreshed independently for each 
mode on average every r = 0.05 L±/v rrns . We define the 
Reynolds number as Re = (Lx/2ir)v rms /p. A typical 
run covers about 50 to 100 crossing times at the largest 
scale. 

As the force renovation time is much shorter than the 
inverse Alfven frequencies of all the excited modes, the 
forcing supplies energies at controlled rates, = \a . 
According to the solution of ||3J| , in this case the + and — 
energy slopes should be independent of the dissipation. 
The results of our numerical simulations are presented 
in Fig. ([1]). They demonstrate that the spectra are 
pinned at the dissipation scale. However, the amplitudes 
of the spectra at large scales are not sensitive to the 
dissipation. As a result, the spectral slopes change with 
the Reynolds number, as to gradually approach the 
balanced spectrum E(k±) oc k^ 2 at large k±. These 
numerical findings agree with the physical expectation 
that large-scale fields are determined solely by the 
large-scale forcing and are independent of the small-scale 



FIG. 2: The field-parallel spectra of the magnetic energy 
(solid line) and the kinetic energy (dashed line) at the field- 
perpendicular wave number k± = 5; Re = 2500, resolution 
512 2 x 256 points. 



dissipation. They however contradict the standard 
model ([2]). In what follows we propose a resolution for 
this inconsistency. 

A model for non-balanced weak MHD turbulence. — 
To derive the model equations, we propose that non- 
balanced MHD turbulence leads to the generation of a 
non-zero average, 



(z+(k) • z-(k')) = 6{k + k')e°{k±)A(k ]l ), 



(5) 



where A(fc||) is concentrated at fey = 0. If the Elsasser 
fields z + and z~ corresponded to Alfven waves, such an 
average would be zero, since Alfvenic fluctuations satisfy 
v = ±b. However, fluctuations at fey = are not waves 
(lu = k\\VA = 0), and the average ([5]) may not vanish. The 
presence of the condensate ([5]) means that the magnetic 
and kinetic energies are not in equipartition at k\\ = 0. 

Physically, non-balanced MHD turbulence is not 
mirror- invariant, as it possesses the non-zero cross- 
helicity, H c — f(v ■ h)d 3 x ^ 0. Non-mirror-invariant 
turbulence can generate large-scale magnetic fields that 
are not in equipartition with the velocity field. This hap- 
pens due to conservation of magnetic helicity in ideal 
MHD, which tends to cascade toward large scales in a 
turbulent state. We propose that such a process is pre- 
served in a driven weak MHD turbulence, although in 
a peculiar fashion - it leads to generation of condensate 
([5]) in the vicinity of k\\ = 0, where the magnetic energy 
exceeds the kinetic energy. Such a condensate is indeed 
observed in our numerics, see Fig. ([2|). 

The dynamics of the condensate are not described by 
the weak turbulence theory, however, they may be ad- 
dressed using certain closure assumptions. It is worth 
pointing out that the same limitation holds for the origi- 
nal system ([2]). As follows from the presence of the (dis- 
function in the integrand, only the e =F ((7|| = 0) modes are 
responsible for the energy transfer. However, if one ap- 
plies Eq. ([2]) to these modes themselves, one encounters 
an inconsistency. The weak turbulence approximation 
is valid when the inverse time of nonlinear interaction is 
much smaller than the wave frequency. The nonlinear in- 
teraction described by the right hand side of @ does not 
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vanish for fe|| = 0, while the linear frequency of the corre- 
sponding Alfven waves, u> — k\\VA, vanishes. Therefore, 
as noted in 0, an additional assumption of smoothness 
of the functions g ± (k^) at kn = was essential for deriv- 
ing the spectra ©■ 

We postpone the discussion of self-consistent conden- 
sate equations for future communications. Here we would 
like to demonstrate the effect that the presence of con- 
densate provides on the turbulence spectra. As we have 
argued, the condensate is expected to alter the cascade 
dynamics in the non-balanced case. Consider the case 
where the imbalance is weak, that is 7 <C 1. In this case 
we expect that the condensate is weak as well. We derive 
the equations for the energies (k) , by proceeding along 
the lines of weak turbulence derivation: we expand the 
MHD equations ([I]) up to the second order in the nonlin- 
ear terms, and split the forth-order correlators into the 
second-order ones according to the Gaussian rule. To 
the first order in e°(k±), the resulting equations have the 
form: 

cV^k) = J M fe , M e^(q) [e^p) - e ± (k)] S( n )d ktPq + 
+A(k\\ ) J Rk, Pq [e k (k±)e°(qx)+e k (q ± )e°(kx)] 4 p9 ,(6) 

where A(fc||) = ReA(fe||), R k , pq = (tt/va)0^± x q_L) 2 (k_L • 
p±)(k_L • q±)/(k 2 _p 2 L q 2 ± ), and d^ pq = 5(h± p± 
1±)d 2 p± d 2 q±. The first term in ([6]) coincides with the 
equation @, while the second term describes the inter- 
action with the condensate. It can be directly verified 
that each of the integrals in ((6|) conserves the Elsasser 
energies E^ = J e ± (k)d 3 k. 

For the stationary solution, each of the integrals in 
((6]) should vanish independently. Equating the first 
integral to zero does not allow one to find the spectra 
uniquely, rather, it leads to the one-parameter family 
of solutions ([2]). Consider the second integral that 
describes the interaction of with the condensate. By 
employing the standard methods of the weak turbulence 
theory, one can demonstrate that the power-law solution 
nullifying the first part of the second integral is unique, 
e°(k±) cx fcj 3 . Analogously, the second part of the 
second integral is zero if e ± (k±) cx fcj 3 . We conclude 
that the presence of the condensate lifts the degeneracy 
of the solutions: the only possible stationary power-law 
spectra of weak MHD turbulence are e^(k±) cx fcj 3 . Al- 
though we do not have the equation for the condensate, 
the above result allows us to predict that in order to 
preserve the scale invariance the condensate should have 
the scaling e°(fc^) cx fc]| 3 . 

Conclusions. — Based on our numerical simulations and 
analytic consideration, we propose that weak MHD tur- 
bulence spontaneously generates a condensate of the 
residual energy E v — Ef, at small fen . We argue that the 



condensate is a consequence of mirror-symmetry break- 
down in non-balanced turbulence. When the turbulence 
is balanced, the energy spectra are E^^ki) cx kj 2 , in 
agreement with the analytic prediction of [J [H, Q . In the 
balanced case the evolution of E^ fields is not affected by 
the condensate. In the non-balanced case the interaction 
with the condensate becomes essential, and we propose 
that no universal power-law spectra exist in an inertial 
interval of limited extent. Both spectra E ± (k±) have the 
large-scale amplitudes fully specified by the external forc- 
ing, and they converge at the dissipation scale. As the 
dissipation scale decreases, the spectral scalings (but not 
necessary amplitudes) approach each other at large k±. 
As a result, the universal scaling kj 2 is recovered for both 
spectra E ± (k±) asymptotically at k± — > 00. 
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